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In this paper, we will develop an adaptive modal identification method for identify- 
ing the frequencies and damping of a flutter mode based on model-reference adaptive 
control (MRAC) and least-squares methods. The least-squares parameter estimation will 
achieve parameter convergence in the presence of persistent excitation whereas the MRAC 
parameter estimation does not guarantee parameter convergence. Two adaptive flutter 
suppression control approaches are developed: one based on MRAC and the other based 
on the least-squares method. The MRAC flutter suppression control is designed as an 
integral part of the parameter estimation where the feedback signal is used to estimate the 
modal information. On the other hand, the separation principle of control and estimation 
is applied to the least-squares method. The least-squares modal identification is used to 
perform parameter estimation. 


I. Introduction 


There has been increasing research interest in wing shaping control technology in recent years. This is 
in response to the need for energy-efficient aircraft design that employs light-weight materials for aircraft 
structures. Reducing airframe operational empty weight (OEW) using advanced composite materials is one 
of the major design considerations for improving energy efficiency. As a result, aircraft wing structures in 
modern transport aircraft have become much more flexible than older-generation aircraft wings. In some 
modern transport aircraft, the wing flexibility can be double that of an older-generation transport aircraft. As 
structural flexibility increases, aeroelastic interactions with aerodynamic forces and moments can potentially 
degrade aerodynamic efficiency and compromise vehicle stability and control. One consequence of increased 
wing flexibility is a reduced flutter margin. Transport aircraft are designed to meet FAA certification for 
flutter clearance. 

Wing shaping control concepts for drag reduction are being studied by NASA to leverage wing flexibil- 
ity for aerodynamic performance.':? By re-twisting a flexible wing and using variable camber aerodynamic 
control surfaces, aircraft wings can have a mission-adaptive capability.? In recognition of the role of aeroelas- 
ticity on aircraft performance and dynamics, NASA Advanced Air Transport Technology (AATT) project is 
conducting research in the area of Performance Adaptive Aeroelastic Wing (PAAW). This research develops 
concepts such as the variable camber continuous trailing edge flap (VCCTEF) to enable wing shaping control 
for aerodynamic performance and dynamics.4 

Structural dynamic interactions with aerodynamics alter natural frequencies and damping of flight struc- 
tures that give rise to flutter.® Flutter is a dynamic stability problem caused by non-conservative work done 
by aerodynamic forces on flight structures. Active flutter suppression control is required to suppress any 
unstable flutter modes that may appear below the flutter clearance speed of an aircraft. Flutter can be 
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modeled by a second-order system with parameter-varying frequencies and damping.® Often times, these 
parameter-varying frequencies and damping are not known accurately due to complex aeroelastic interac- 
tions. Real-time modal identification can potentially helps to determine the unknown parameters of the 
aeroelastic model of an aircraft. The information obtained from the modal identification can be used to 
devise an active flutter control to improve the flutter margin. 

In this paper, we will develop an adaptive modal identification method for identifying the frequencies and 
damping of a flutter mode based on model-reference adaptive control (MRAC) and least-squares methods. 
The least-squares parameter estimation will achieve parameter convergence in the presence of persistent 
excitation whereas the MRAC parameter estimation does not guarantee parameter convergence.” The sep- 
aration principle of control and estimation can be applied to the least-squares method. The least-squares 
modal identification is used to perform parameter estimation. Then, an adaptive flutter suppression control 
could be designed using the parameter estimates obtained from the least-squares modal identification. The 
least-squares parameter estimation method is demonstrated on an airfoil flutter problem. The recursive 
least-squares approach is observed to be highly effective in achieving parameter convergence. 


II. Flutter Model 
The aeroelastic equation that describes flutter of an aircraft wing is given by 
Mg+Cq+Kq=Q (1) 


where q € R” is the generalized coordinate vector, M € R”*” is the generalized mass matrix, C € R"*” is 
the generalized damping matrix, K € R”*” is the generalized stiffness matrix, and Q € R” is the generalized 
aerodynamic force vector. 

The generalized aerodynamic force vector Q is a complex function of the aerodynamic characteristics, 
the so-called reduced frequency, the generalized displacement q, the generalized velocity q, the generalized 
acceleration g, the control surface deflection 6 € R™ which comprises m control surfaces, the deflection rate 
é, and the deflection acceleration 6. So, in general, it can be expressed as 


Q= Qa (Poo, Voo) q+ Q4 (0355 Voo; k) q+ Qa (Poo, Voo, k) q+ Q5 (pissy Voa)0 
+Q5 (Poo, Voo, k) 6 + Qs (Poo; Voor k)5 (2) 


where k = - is the reduced frequency, w is the flutter frequency, ¢ is the mean aerodynamic chord, V3, is 


the airspeed, and py is the air density which is a function of altitude. 
The generalized coordinate of a flutter mode can be described by qf (t) = qo; e(—Svtiw)t where ¢ is the 


modal damping. Flutter occurs when the damping of a flutter mode crosses zero from ¢ > 0 to ¢ < 0. 
T 
Let q = af 4s where gy € R is the generalized coordinate of the flutter mode and gq, € R"~' is the 





generalized coordinate vector of the remaining modes. 
The aeroelastic equation for the flutter mode can be separated from Eq. (1) as 


(My — Qya,) d¢ + (Ce — Qra,) d¢ + (Kp — Qa) a¢ = (Qra, — Mye) Ge + (Opa. — Cjs) ds 
+ (Qyq, — Kgs) ds + @ p39 + Q p55 + Qp56 (3) 


(M, — Sean) Os ae (Cs = sda). Ge = (Ks = Qsq.) Us = (Qsas a M,1) af ate (Osa) ~~ Cs) af 

+ (Qsa, — Kap) Up + Q,39 +. Q,59 + Qs55 (4) 
Ky—-Qpra; Cr—-Qha;z 
: Mys—-Q sa; Ms—-Qra;z” 
6 = doe’ where w = wy + € is the excitation frequency with € as a small number. The response of the 
generalized coordinate of the flutter mode is described by q¢ = qo; e(“t-$) where ¢ is a phase lag, while the 


Let we = and 2¢pwy = Suppose the control surface deflection is commanded by 


response of the generalized coordinates of the remaining modes is described by gs = qo, et) where 7 is 
the phase lag. Then the responseq, is obtained as 


ds = F-'Gas + H's (5) 
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where F = —(M, — Ssg,) w? + (Cs — Qsq,) iw + (Ks — Qsq,), G = — (Qea; — Ms) w? + (Quay _ Css) tw + 


(Ques _ Kat)s and H = —Q,3u? + Q giv + Qs6. 
The response of the flutter mode is computed to be 


[— (Qpa. — Myps) w? + (Qpq, — Crs) iw + (Qpq, — Kps)] (F~'Gar + H~'6) 
(My — Qya,) (-w? + 2¢ purpiw + w?) 





df = 


(-2 jg” + Qygiw + Qz6) 6 


(My — Qsa;) (-w + 2G pw pi + w?) 


Since w = wy +e, ¢ + 0 ase — 0. Therefore, the amplitude of gy is larger than that of q,. Thus, the 
flutter mode is the most dominant response in the presence of persistent excitation input from 6. Therefore, 
it can be expressed as . ; 

Gp + 2 pwr de + w5qe = b16 + bod + b36 + Ag (7) 
where A, is an unmodeled dynamic due to the remaining modes. 


Suppose the measurement is an acceleration of the wing bending deflection at a given location z. Then, 
the acceleration is expressed in terms of the generalized coordinates as 


w= O(z)G= of (2) Gr + ®s (y) Gs (8) 


where ®(z) = br (z) ,(z) | € R” is the mode shape vector, ¢f (z) € R is the mode shape at the 


measurement location for the flutter mode, and ®, (z) € R”~! is the mode shape vector for the remanning 
modes. 

When the excitation frequency approaches the flutter frequency, the acceleration response is predomi- 
nantly due to gf. Thus the acceleration can be written as 


w = bf (z) dp + Aw (9) 


where A,, is the contribution of the remaining modes to the acceleration measurements which is assumed to 
be small. 
Therefore, the measurement equation is approximated as 


W + Wpwph + ww = c1d + cod +036 + Ag (10) 


where c; = ¢y (z) bj, 1 = 1, 2,3. 

We will use this equation to estimate the fluter frequency and damping. Then, using this information, 
an adaptive control can be designed to suppress the flutter mode. 

Example: Consider an airfoil in pitch and plunge as shown in Fig. 1. 





Figure 1. Airfoil in Pitch and Plunge 
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Consider a quasi-static flutter problem whereby we neglect the aerodynamic mass and damping matrices. 
The equations of motion are written as 


La] [ sass Aas || | 


m Meg Vv : 
0 269 ko — AY Stadenc 








MEcg I 6 6 
G(k) C1, qc00” 
4 | kn F (k) c1,,doo€ A} _ —F (k) c15qooe 5 SS sls 
0 ke —F(k)ci.doce || 0 F (k) (ems + £¢1s) doo G(R) (ema tet, Joos” 


Given m = 1.86733 slug/ft, I = 15.5611 slug-ft, kj, = 500 Ib/ft?, ke = 10,000 lb, Ch = Cp = 0.01, c= 10 
foe U15G eg = Oe = 7 2m __ for subsonic flow, Mj. = 0.4, and poo = 2.37756 x 107° slug/ft? at 








1—M2, 
standard sea level. The control derivatives are given by 
Ci, = 2 cos (a-2%) +2,/2 (1- “)| (11) 
T c c c 
cms =~ $2 [eos (1-24) +22) JL (1-2) (12) 
27 c c c Cc 


where cr = 2 ft is the flap chord. 

Figures 2 and 3 show the frequencies and damping of the pitch and plunge as a function of the equivalent 
airspeed in knots (KEAS). The flutter mode corresponds to the plunge mode which occurs at an airspeed of 
61.7353 knots and a frequency of 16.2804 rad/sec. The damping plot actually shows the structural damping 
required to achieve neutral stability. It is two times the negative value of the viscous damping ratio. 


























30 T T T T I I 
—— Plunge Mode 
—— Pitch Mode 
Flutter Point 
8 257 
ToT 
& 
3 
=] 
g 
oa 
id 
Ww 20+ 
15 1 i 1 i 1 i 
0 20 40 60 80 100 120 140 


Airspeed, KEAS 


Figure 2. Frequencies of Airfoil in Pitch and Plunge 
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Figure 3. Damping of Airfoil in Pitch and Plunge 


III. Parameter Estimation of Second-Order System 
Consider a second-order system with matched uncertainty in a state-space form 
é=Ar+Blu+0*'&(2)| (13) 


where A and B unknown but sign of b is known. 
Assuming that there exist K, and k, that satisfy the model matching conditions, and furthermore that 
Am and B,, also correspond to a second-order system. Then A and B can be identified. Let 











a-|°, 2 |.2-[° el 1 fe-| 2 a 
—wi —2Cw, b —we, —2mwm Dm 
The model matching conditions are 
A(t) + B(t) Ky (t) = Am (15) 
B (t) ky (t) = Bm (16) 


Ke= (8°) "8 (4m A)= 5 [0 6] ] soe au an | 





= : =W2 Ge —2C aus, +26; |) 1) 
c= (8°8) BB = [0 #]|,) |= % a 


where A, B , Wn, and ¢ are estimates of A, B, wn, and C¢, respectively. 


Let A(t) = A(t)— A and B(t) = B(t)—B be the parameter estimation errors. Now rewriting the plant 
model as 


ie (4-4) a4 (6-8) [u+ 0°78 (x) (19) 
The adaptive control is designed as 
u=K,(t)a+k,(t)r—O!' (t) ®(z) (20) 
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Then, substituting Eqs. (15), (40), and (17) into Eq. (18) yields 
r= (4 = A) c+ B [Kya + ker — @'® (x) + 0*' © (x)] — B [Kya + ker — O18 (x) + O°7 8 (2)| 


= | A+BK,—-A| «+ Bk,r—BO'®(r)— B(Kyt+kpr) (21) 
—_——~—— “ 


Am Bm 


where © (t) = © (t) — ©* is the matched uncertain parameter estimation error. 
Then, the tracking error equation is established as 


6=tm—&£=Anet+Ar+ But BO &(z) (22) 


where u = K, (t)x +k, (t)r. 
Choose a Lyapunov candidate function 


bBo 2 2. B'B 43 ~ 
4 (c. A, B, 6) =e! Pe + trace (473747) ++ Fy. + |b} O'TS'O (23) 
b 


where [4 =) > 0 andTe = i > 0 are positive definite adaptation rate matrices. 
Vy (c. A, B, 5) is evaluated as 


. ~~ ~ ~ ~ ~T 2Bt B im 
V (c,4, 8,6) Se Osteet Pp [Av + Ba + BO ®(c)| 4+ trace (2413747) + = td 67rg16 (24) 


The following parameter estimation adaptive laws are then obtained 


A= —Tyze'P (25) 
B = —yPeu (26) 
0 = —l'e® (x) e' Psgn (0) (27) 


where P = £4. 
It follows that e (t), A(t), B(t), and © (t) are bounded since 


V (c,,B,6) = ~e7Qe < —Amin (Q) lel? (28) 


V (c. A, B, 6) has a finite limit as t > oo since 
V (t 3 cw) = V (to) -f[ Amin (Q) |lel|” dt < co (29) 
to 


It can be shown that V (e, A,B te) is uniformly continuous because V («, A, B,®) is bounded. Then, 


applying the Barbalat’s lemma, one can conclude that the tracking error is asymptotically stable with 
e(t) + 0 as t > oo. However, this does not imply that the parameter estimation errors A (t), B(t), © (t) 
tend to zero. So parameter convergence is not guaranteed with MRAC. 


Let 


A=[0 1 |A= (30) 


in 
— 
Oo 
Hm 
— 
rr 2 
| 
> oO 
wo 
| 
bo 
I™> 
& 
3 
[_ | 
lI 
— 
- 
3nNu 
| 
bo 
ms 
E€ 
3 


and since 
b=|0 1|B=[0 ar) (31) 
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then the adaptive laws can be expressed in terms of the estimates of the unknown quantities w,,, ¢, and b as 














A’ =-Tyze'™P cae —Tyze'P (32) 
ie — 0 1 Peu = —»P! eti _ —ypiie' P (33) 
Let 
ra=|% 9 | s0 (34) 
YC 
Then d 
4 (a2) = yore" P (35) 
or a 
tel P 
oe (36) 
d ro ~ 3 et Tp 
Hi (—2¢an ) = — Qing — 2¢dn = —Yyet2e P (37) 
or 


(ncaatn _ wid) e'P 


= 38 
( ms (38) 
To prevent the possibility of w,, = 0 or b = 0 that will cause the adaptive laws to blow up, both the 


adaptive laws for estimating w,, and b need to be modified by the projection method as follows: 





qwtse Pigg 0 or if @, = dé» >0 
oie Bo, I Wy ie or 1I Wy Wo and Wyn = (39) 
0 otherwise 
A —yie'P if 6 > bo or if 3 = bp and dt} >0 
b= (40) 
0 otherwise 


In the modified adaptive law for w,, it is assumed that w, is always a positive quantity for a physically 
realizable system. 
Example: For the following second-order SISO system 


#1 + Cw, + wre, =b [u +0*'® (x)] 


where ¢, w,, and b are unknown, but b > 0 is known, we design an indirect adaptive controller using the 
following information: ¢, = 0.5, Wm = 2, bm = 4, r (t) = sin 2t, and 


For simulation purpose, the unknown parameters are given to be € = —0.5, w, = 1, b = 1, and O*! = 
0.5 —0.1 , and all initial conditions are assumed to be zero, except for w,(0) = 0.8 and b(0) = 0.6. For 
simplicity, we use the unmodified adaptive laws for w, and b. Use Ywo = Ye = Y = 10 and Te = 10/. Figure 
5 shows the closed-loop response of the second-order system to MRAC. All signals are bounded and the 
tracking error tends to zero asymptotically. Figure 5 shows the parameter convergence which is quite good 
even though the excitation frequency w = 2 rad/sec is not too close to the plant frequency w,, = 1 rad/sec. 
The adaptive control is able to stabilize the unstable plant. 
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Figure 4. Response to MRAC of Second-Order System 
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Figure 5. Response to MRAC of Second-Order System 


Note that the adaptive modal identification based on MRAC assumes zero contribution from the control 
surface deflection dynamics due to the terms 6 and 6. This may be a reasonable assumption if the excitation 
frequency is small. Increasing the excitation frequency can cause these terms to become significant. In such 
a case, the dynamics of the control surface deflection must be accounted for. Moreover, we will consider 
least-squares parameter estimation method which is generally more robust than MRAC and can guarantee 
parameter convergence if the persistent excitation condition is satisfied. 

Going back to the measurement equation (10), we establish an estimator model as 


tha + 2a pb + O2w = G10 + God + 635 (41) 


The output is the acceleration w. We compute the plant modeling error between the true output and 
the estimated output wg as follows: 


€= tig — tH = 2H pwd — Dw t+ Gb+ d+ GS (42) 


Note that the system is open-loop with the control surface deflection providing the persistent excitation 
signal 6 = doe*“?. The parameter estimation will be computed by a least-squares gradient method by 
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minimizing the cost function J = se. Taking the gradients of the cost function with respect to the 


parameter estimates yields 














SION ot ie (43) 
a (a) 
es = —we (44) 
(2) 
OJ + 
oe 4 
Oe, de (45) 
OJ: 
=4d 46 
dé m6) 
OJ 
=4d 47 
dé; oe 
The least-squares gradient adaptive laws for modal identification are then obtained as 
d xO OJ 
dt (oF) 9 (a) “Yur WE ( 
w@ 
f 
d/s. OJ : 
a (61) = ee = Dre (49) 
0 (as) 
* OJ es 
c= “See Rs = —Yc, 0€ (50) 
2 OJ : 
2 Tea Be = =Vex0€ (51) 
: OJ 
C3 = TS Bes = —Yc, 0€ (52) 


where y,, > 0, yc > 0, and 7, > 0,7 = 1,2,3 are the adaptation rates. 
The frequency estimation adaptive law can also be written as 





2 WoW 
= 53 
Eon, (53) 
2 Ayew we — yCwe 
C= oe (54) 


In the presence of the persistent excitation conditions + i wdt = Il) > 0 and af 67dt = (15113 > 0, 


2 2 ‘, 
exponential parameter convergence will be achieved. Since 6 = dge**, then \|51)3 = 0: and w = woet-9) , 


2 
then I|wl|5 = <2. Therefore, the least-squares adaptive modal identification will achieve parameter con- 


vergence. Then, in theory, this implies &; (t) > wy, ¢(t) > Cy, &1 (t) > c1, 2 (t) > ca, and é3 (t) > cs 
in the limit as t — oo. However, in practice, the parameter convergence is only approximate because the 
estimator model does not represent the true plant. The modal contribution of those modes not captured in 
the estimator model will influence the parameter convergence. 

Adaptive flutter suppression control can be designed based on the parameter estimation. We assume a 
time-varying feedback law 6 (t) = ky (t) w (t) +r (t), where ky (t) is an adaptive feedback gain and r (t) is a 
command signal which is used to reject time-varying disturbances. The adaptive flutter suppression control 
is to achieve a stable reference model that represents the closed-loop dynamics as follows: 


Wm + mmm + w2,wm =0 (55) 
while the real plant is represented by 
th + Wwyd + wew = ¢16 + cod + c36 + v(t) (56) 
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where v (t) is a time-varying disturbance. 
We reply on the separation principle of control and estimation to use the estimator model as if it is the 
true model for designing a controller. Then, upon substitution, we get 


1 — &ky) t+ (26 pd" — Werky — Coky ) w+ (02 — Gky — Cokw — esky) w=cartettertu (57 
ft f 


Thus, the adaptive law for the feedback gain k,, is obtained as 





Re 2C pp — Win + 2Xmwmerkw ~ Cakuw (58) 
2¢1 
¢mWm is chosen to be such that 2¢mWmé1 — €2 < 0 if @) > 0 or 2¢nwmer — C2 > 0 if G <0. Then, k,, (¢) 
tends to : 
Tr 2 mem ~~ 2 A 
frome Aes (59) 
2CmWme1 — C2 
in the limit as t > oo. 
This implies the damping ratio tends to 


Zee (60) 





The command signal to cancel out the disturbance can be designed by constructing a disturbance esti- 
mator as 


b = ~ Te 5 = VE (61) 
Then, the steady state command for disturbance cancellation is computed as 


c3r = 0 (62) 


Thus, the adaptive command signal could also be computed from the least-squares parameter estimation 
as 
f= —r€ (63) 


Example: Consider the previous example of an airfoil in pitch and plunge. We will consider an open-loop 
modal identification of the plunge mode by applying a swept sine excitation signal from 15 - 17 rad/sec. 


6= 0 sin (wy = wo) = + Wo t (64) 
? 


The plant modeling error € is computed from the estimator model and the accelerometer output located 
on the airfoil. The top of figure 6 shows the plot of the estimated frequency of the plunge mode. The 
adaptation rates are chosen to be y%, = 200 and ye = 1. The plant corresponds to a natural frequency 
of wr = 15.908 rad/sec and a damping ratio of ¢ = 0.0036. It can be seen that the estimated frequency 
approximates the true frequency of the plunge mode very well. 

The lower plot of figure 6 is shows the estimated damping. The signal is noisy as expected because the 
true value is a very small number. Nonetheless, one can compute a mean value of about 0.0035 between 12 
sec and 20 sec. This also agrees well with the true value. 

Figure 7 is the plot of the estimated output from the estimator model and the true output. The estimated 
output converges to the true output after about 8 sec. 
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Figure 6. Estimated Frequency and Damping of Plunge Mode with Least-Squares. Swept sine wave from 
15-17 rad/sec. 
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Figure 7. Estimated Output of Plunge Mode with Least-Squares 


One difficulty with this approach is the sensitivity to the excitation signal. For example, if less is known 
about the parameter, and a broader sine excitation signal that varies in frequency from 10 - 20 rad/sec is 
used, the parameter estimates can wander away from their true values as shown in figure 8 below. 
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Figure 8. Estimated Frequency of Plunge Mode with Least-Squares. Swept sine wave from 10-20 rad/sec. 


Next, in an attempt to compactly excite a broader range of frequencies, a 2 Hz square wave is used as 
the excitation signal 6. For é,a signed Kronecker delta is used at the square wave axis crossings. jy, is set to 
10. Lower frequency square waves (below 10 Hz) produced the best results, and higher amplitudes decreased 
convergence time more stably than increasing y,, as shown in figure 9 below. 
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Figure 9. Estimated Frequency of Plunge Mode with Least-Squares. Various square wave amplitudes at 2 Hz 


FlIgure 10 shows that one downside to this method is that the estimated plunge mode damping is even 
noisier than before. 
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Figure 10. Estimated Damping of Plunge Mode with Least-Squares. 5 degree square wave at 2 Hz 


Finally, random Gaussian noise is used to excite this Least-Squares adaptive system. Using adaption 
rates of y,, = 1000 and ye = 1, and a 10 Hz signal with zero mean and one degree standard deviation, both 
paramerts are shown to converge within the neighborhood of their true values in figure 11 below. However 
both remain noisy for obvious reasons. 


18 


Estimated Frequency of 
Plunge Mode [rad/sec] 





8 1 1 1 1 1 1 1 1 1 J 
0 5 10 15 20 25 30 35 40 45 50 


Time [sec] 





¢ 
Oo 
= 





Estimated Damping 
of Plunge Mode 
oO 


0 5 10 15 20 25 30 35 40 45 50 
Time [sec] 


So 
o 
© 


Figure 11. Estimated Frequency and Damping of Plunge Mode with Least-Squares. Gaussian noise with zero 
mean, one degree of standard deviation at 2 Hz 


Because of the signal-to-noise ratio for the parameter estimation of the damping ratio, the recursive 
least-squares method is used where the adaptation rate y,, and ye are modified as follows: 


Ju = —Nuyw (65) 


Vo = —neyew? (66) 


where 7, > 0 and 7¢ > 0 are the tuning parameters, chosen to be 0.5. 
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The plots in figure 12 show the estimated frequency and damping of the plunge mode generated by the 
Recursive Least-Squares method using the same excitation signal and initial adaption rates of those used 
in figure 6. The estimated frequency converges to a value of 15.933 rad/sec versus the true value of 15.908 
rad/sec. This is an excellent agreement. The damping value also trends toward convergence with the final 
mean value of 0.0044 versus the true value of 0.0036. Compared to the results of figure 6, both the estimated 
frequency and damping stay within a tighter bound and exhibit less noise. 
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Figure 12. Estimated Frequency and Damping of Plunge Mode with Recursive Least-Squares 


IV. Conclusions 


This paper presents a modal identification approach to estimate the frequency and damping of a flutter 
mode of an aeroelastic system. The parameter estimation is based on the least-squares method using the 
output error. An application for an airfoil flutter is studied. The method shows that the parameter estimation 
is highly effective in modal identification of an aeroelastic system. In particular, the recursive least-squares 
demonstrates to be more effective than the standard least-squares gradient method in suppressing noise for 
signals with low signal-to-noise ratio. Using the adaptive parameter estimates, a flutter suppression can 
be designed to stabilized a flutter mode using the separation principle of control and estimation. In future 
work, we will develop flutter suppression control based on the modal identification method as presented in 
this paper. 
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